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a formulation involving three auxiliary composite fields. The effective potential is constructed at 
leading order in the 1/n expansion for a condensate belonging to the center of the U(n) group. A 
wide region is found in the coupling space where in addition to the condensate proportional to the 
unit matrix, also metastable minima exist even in the chiral limit, in which a further breakdown 
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conjugate to this component of the order parameter changes this state into the true ground state of 
the system. 
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The numerical details of the pattern of chiral symmetry breaking in theories of three-flavored strong interactions 
are essentially dictated by the observed mesonic spectra PHI] . Potentials of effective models are constructed with the 
requirement that, in the ground state, a condensate proportional to the 3x3 unit matrix should arise spontaneously 
and split the degeneracy of the parity partner states. This fundamental structure is slightly modified by an external 
Q_i source also proportional to the unit matrix, which generates mass for the pseudo-Goldstone pseudoscalar mesons. The 
presence of a further external field proportional to the Tg generator results in the kaon-pion mass splitting via the 
r"" 1 1 formation of a supplementary condensate along the same direction of the Lie algebra. Finally, the effect of the axial 
anomaly is invoked for interpreting the emergence of a heavy 77 meson. A complete parametrization of the three-flavored 
CN| linear sigma model has to rely also on some mass information concerning the scalar sector. 

It can be expected that this pattern will not remain unchanged when the number of flavors is increased. The 
investigation of such a sequence of theories is employed when looking for a reliable starting point in constructing the 
nonperturbative solution of strongly coupled effective meson theories. The solution of the large n extension of the 
^-j- two-flavor theory from an 0(4) symmetric theory over an O(n) symmetric one is well-known [5] but does not bring 
• any new structure in the condensate since only a unique type of vacuum condensate might emerge [up to global 0(n) 
^Cj rotations]. Anyhow, this extension found several applications in modeling chiral symmetry breaking in QCD [7j 
I and also in technicolor theories 019]. Applications to nonequilibrium phenomena, e.g., to the quantum dynamics of 
disoriented chiral condensates [TUHH], also benefit from this extension, 
r* More complex symmetry breaking patterns might arise when the three-flavor U (3) x t/(3) symmetric effective meson 
. £^ theory is generalized to large n. The implementation of the restricted large n extension described above retains only 
the 0(2n 2 — 18) invariant part in the quartic potential. Here, the main obstacle in front of wider applications is the 
I— 1 missing leading order solution when n — > 00, which would make explicit the differences between the ground states of 
^ the U(n) x U(n) symmetric and the 0{2n 2 ) symmetric models. A first step in this direction was made with the matrix 
generalization of the Hubbard-Stratonovich transformation involving the complete quartic potential of the U(n) x U(n) 
symmetric model [13] . Based on this rewriting of the model a saddle-point solution was constructed for the field theory 
defined with the auxiliary matrix field |I4| . It was conjectured that this solution describes the leading order behavior 
of its ground state for a large flavor number. The approximate solution was applied directly to the thermodynamics of 
the three-flavor effective meson model (see also [15] ). In the actual solution of the matrix saddle-point equation, the 
authors restricted their Ansatz to a condensate proportional to the unit matrix, which reduced once again the quartic 
potential to its 0(2n 2 ) invariant part. As a result the part of the vacuum condensate pointing in the algebra space 
along the Ts direction remained proportional to the T% component of the explicit symmetry breaking external source. 
To our actual knowledge no published further step was attempted towards the exact large n analysis of the U(n) x U(n) 
symmetric matrix model until the proposition |I6] appeared. 

The new feature of |16j was the introduction of two auxiliary fields: beyond an 0(2n 2 ) invariant which corresponds 
to earlier Ansatze, also a U (n) vector appeared with nonzero saddle-point value. Following the introduction of these 
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auxiliary fields, the large n limits of the one-point and two-point equations were found and solved with a physically 
arguable supplementary restriction. Although only a condensate proportional to the unit matrix (v = vo x I) was 
assumed, the introduction of the new composite field with nonzero saddle value led to new features not present in 
the large n solution of the 0(2n 2 ) model (for instance, nondegenerate spectra emerged in the scalar meson sector). 
The supplementary assumption invoked self-consistently during the construction was the heaviness of the scalar fields 
relative to the pseudoscalar degrees of freedom. A large domain was found in the three-dimensional coupling space 
formed by the two independent quartic couplings and the external source h = ho x I, where this self-consistency is 
fulfilled. 

The present paper is a continuation of [16] and expands the discussion of the proposed solution in three new directions, 
shortly described below. 

A natural next step towards the phenomenological application of this solution is the extension of the condensate by a 
term proportional to the "longest" diagonal generator: T^ lag ~ diag(l, 1, 1, — n+1), which coincides, for n = 3, with 
Tg. There is the possibility for introducing a more general condensate belonging to the Cartan subalgebra of U(n), but 
the more restricted Ansatz T Q v + T^ a9 vf a9 is the one smoothly connected to the relevant n — 3 case. In the present 
study, we construct an approximate large n effective potential at leading order for such a condensate and search for its 
global minimum with and without the presence of the corresponding explicit symmetry breaking source, T^ ia9 h^ a9 ' . 

The second new feature is the replacement of the single composite U (n) vector introduced in [16] by the sum of two 
fields, each of definite parity. One of them mixes with the scalar, the other with the pseudoscalar elementary fields. The 
composite pseudoscalar is assumed to be similarly light, like the elementary 7r a ; therefore, its dynamics will be included 
into our analysis, while, in [T5], the sum of the two fields dominated by the scalar composite was classified to be "heavy," 
and its fluctuations were omitted from the contributions to the equation of state and the saddle-point equations. In 
the present refined analysis, the mixing in the pseudoscalar sector results in a spectrum, where also heavier excitations 
("kaons") might appear naturally. The self-consistency of this solution requires that even the heavier pseudoscalars 
should be considerably lighter than the lightest scalar excitation. We call the part of the coupling space where this 
condition is fulfilled the "allowed space." 

It will be shown that, in the allowed coupling space, besides the "trivial" solution (e.g., vq ^ 0, v^ la9 = 0), also a 
nontrivial local minimum of the effective potential with w ^ 0, w* aff ^ exists even for h = h^ ag = 0. We shall 
explore the behavior of this minimum also for ho ^ and find that there is an upper bounding value, beyond which 
the nontrivial local minimum disappears. 

In order to decide which is the global minimum (i.e., the ground state), we shall map out the effective potential along 
the v^ a9 axis. The construction of this potential and the proof of its renormalizability is the third new result of the 
present investigation. First, we interpret the behavior of the potential for ho ^ 0, h°^ a9 = 0. We find that the trivial 
solution is favored energetically everywhere in this subspace. In this way, we demonstrate the global stability of the 
solution presented in [16] in the absence of explicit symmetry breaking along the T*" 9 direction. In some points of the 
coupling space, where the difference in the potential energy compared to the trivial minimum is relatively small, also 
the effect of a nonzero h^ a9 is explored. Dialing this external source one can tune the state with v°l ia9 ^ to become 
the true ground state of the model. 

In Sec. II, we shortly reformulate the model with three auxiliary fields. In Sec. Ill the renormalized spectra of 
the pseudoscalar sector are worked out, since in view of the assumed mass hierarchy, only the fluctuations of these 
fields contribute to the equations of state and the saddle-point equations. The effective potential is renormalized and 
evaluated in Sec. IV with the pseudoscalar propagators determined in Sec. III. The numerical solution covering the 
part of the coupling space where the assumed mass hierarchy is self-consistently valid is presented and discussed in Sec. 
V. Our conclusions will be summarized in Sec. VI with some hints towards the phenomenological consequences of the 
new type of ground state described in the paper. 



II. FORMULATION OF THE MODEL WITH AUXILIARY FIELDS 



The Lagrangian of the model can be written with the help of n 2 pseudoscalar (7r a ) and n 2 scalar (s°) components as 

L = \[{d^ a ) 2 + {d^ a ) 2 ~m 2 {{s a ) 2 + {K a ) 2 )]+^h a s a 

- £^((^ a ) 2 + (^) 2 ) 2 -£[(^ a ) 2 + (^ a ) 2 ], a) 

where the U(n) vectors i = 1, 2, arc defined as 

ttcl jabct b c , __b„c\ jja rabc be fc%\ 

Ul = -d [S S + 7T 7T ), U 2 = } S 7T , (2) 

with d abc and f abc being the symmetric and antisymmetric structure constants of the U(n) group, respectively. We 
note that, in [16], the sum U1 + U2 was treated as a unique combination, and only two auxiliary fields were introduced. 
In view of U^U^ = 0, the square of the sum becomes the sum of the two squares as it has already been applied above. 
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Correspondingly three auxiliary composite fields are introduced by adding the following constraints to the Lagrangian: 

2 , / / \ 2 , , , x 2 



AL = -~| X-i 



(3) 



In the sum C = L + AL, a symmetry breaking pattern corresponding to the shifts 

s a s a + V2^v q d qa , US US + nV2v q f aqc n c , 
US -* US + nV2d abq s b v q + n 2 d apq v p v q . 



(4) 



is realized. The indices "q" denote the diagonal generators of U(n). The index q = refers to the unit matrix, and 
q = 2, n refer to the elements of the Cartan subalgebra of SU(n). 

After the introduction of the auxiliary field variables and the shifts, the full Lagrangian has the following form: 



L = 



(fl„s«) + (d M vr a ) 2 - m 2 [2n 2 v 2 q + 2V2n 2 v q s« + {s a f + (7r a ) 2 Jj - - (X 2 + (Yff + (Y 2 a ) 2 ) 
/ 2m 2 h q {s q + V2n?v q ) + i^^X (2n 2 v 2 q + 2V2n?v q s q + {s a ) 2 + (^ Q ) 2 ) 



+ i 



Y?U? + nV2v q d abq Y 1 a s b + Y?n 2 d apq v q v p 



Y 2 a US + nV2v q f aqc Y 2 a Tr c 



(5) 



The background v q introduces mixing of the pair s a ,Y" for every value of the index a. In addition, it induces nonzero 
values for YJ , and also a (2n + l)-dimensional mixing sector is being formed for the propagators of these degrees of 
freedom: s q , X, and Y q . On the other hand, pions will mix (for some index values at least) with the composite field Y 2 a , 
which has negative parity. Since the ground state of strong interactions respects parity, we expect to have a quantum 
solution of the model, where Y 2 = 0. Note that, if Y 2 had nonzero value, scalar and pseudoscalar fields would mix, 
and the analysis of the mass spectrum would be a lot more complicated. Instead of going through the most general 
analysis, we assume from the start that 7r° and s a fields do not mix and, below, we show that this is consistent with 
1? = 0. 

We construct correspondingly a quantum solution, where the saddle-point values of X, Y q are nonzero, Y 2 = 0, and 
the two-point functions corresponding to the above-mentioned coupled (mixing) sectors are assumed not to vanish. In 
our approximation, only the quantum fluctuations of the light 7r-fields enter the equations of state of the condensate 
and the saddle-point equations; therefore, we fully avoid, in this paper the analysis of the effect of the more complex 
symmetry breaking pattern on the s sector. The tt sector breaks up into two separate groups of fields. Those for which 
fuqw 7^ form two-dimensional subsectors with the corresponding fields [Yj",^]. Their masses depend on the actual 
background configuration explicitly. The masses of the pions which do not mix with Y 2 (because f uqw = 0) depend 
only on the saddle-point values of the auxiliary fields. 

In the following, we focus our attention on the case where the condensate consists of two components, e.g., q = {0, n}. 
Even though this paper is about a strict large n study and, therefore, no specific reference for the n — 3 case will be 



given, for the sake of brevity, the notations T^ lag = Tg, v°\ 



diag 



f'8 



, and hfS ag — h$ will be used. 



III. EQUATIONS OF STATE, SADDLE-POINT EQUATIONS, AND PSEUDOSCALAR PROPAGATORS 

Standard rules [171 of constructing the derivatives of the quantum effective action T with respect to the fields can be 
readily applied to km . The contribution involving propagators of the heavy scalars are simply omitted in view of the 
assumed mass hierarchy. 

In the equations of state (EoS) and the saddle-point equations (SPEs), the useful combination M 2 can be introduced: 



M 2 = m 2 - -{^fh^X + v/^Yj ). 



(6) 



The SPEs of Y x , Yj, and X and the EoS of vq and v$ take the following form: 





.^252 

= i 




n 


Y? 


i f 9 ^ 




Y n 





= nV2(h 





= nV2(h 



n 2 (v 2 +V 2 8 ) + ^ J G n a n a(p) 

1 



2nV2nv v 8 + n d 8SS v$v s + -d Sab J G„ a ^b(p) 



r2g 2 X = y^Y?, 



nV2(h - M 2 v ) + i2^Y*v 8 



1Kb 



GW&(p). 



(7) 
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Here, propagators G^ a ^b and G^ay* refer to the elements of the complete 2n 2 x 2n 2 propagator matrix in the tt — Y% 
sector. The structure of this matrix (e.g., the mixing of the fields) will be made explicit later in this section. From the 
above equations, it is clear that the mass splitting arising in the pion sector plays an important role when one evaluates 
the tadpole-type contributions. As for the saddle-point equation of the auxiliary field Y2, we obtain Y 2 = 0, since 
the only contributing quantum fluctuations would come from the tadpole of the scalar-pseudoscalar mixing propagator 
G s a^b, which was previously excluded from the solution on the basis of the parity conservation of strong interactions. 
This shows that the F 2 a = assumption is consistent. 

For the various two-point functions of the tt sector, it is convenient to introduce the following inverse propagator 
commonly appearing in several two-point functions: 



iD^(p)=p 2 



M 



f2 



(8) 



We substitute for all three-point functions their classical expressions, which can be read off from ^ (and listed explicitly 
in [16J), and close the coupled set of Dyson-Schwinger equations at the two-point level (Bare Vertex Approximation, 
i.e. BVA). However, it should be clear that this solution nontrivially differs from the application of BVA to the original 
formulation without auxiliary fields. 

First, we discuss the potential mixing of the ir a fields among themselves. Using BVA at the level of the auxiliary 
field-meson-meson vertices and omitting the bubble contribution fully, since at least one heavy scalar line is involved, 
one finds for the inverse propagator matrix element of two pion fields 



iD n 6 WU + i 



(9) 



The only mixing connects the components and 8, which scales, however, for n — !• 00 as 1/n and is, therefore, 
neglected in what follows. The mass of all diagonal components stays with M 2 , except the field "8": 



M. 



88 



(10) 



A large set of tt fields that consists of coefficients of the nondiagonal generators might mix with the composite fields 



Y 2 , since 



T, 



-i^2ng 2 f u8 



(ii) 



The structure constants f uSw are nonzero only for off-diagonal T U \T W generators when [T U ,T W ] is diagonal. Since, 
in the off-diagonal class, the generators are proportional to the "transmutations" of the Pauli matrices a x ,<r y with 
nonzero matrix elements shifted to the row-column positions (j,k), one might refer to them using the extended index 
notation (x, j, k), [y, j, k). One quickly recognizes that each of the 2(n — 1) pairs u — (x,j < n,n),w — (y,j < n, n) and 



(y, j < n, n), w = (x,j < n, n) form coupled 2x2 sectors of tt u , Y™ . In the limit 



one has /( 



— l/v2 + C(l/n). Below, wherever tt fields with these index values u, w occur, we simply use the index K with reference 
to kaons. Also note that, starting from the definition of the structure constants, one finds that the only nonzero d 
coefficients at n = 00 involving off-diagonal generators are rf(a;j,n),(Kj',n),8 = ^(y,i,n),(j/,j,n),8 

= -l/ v / 2 + 0(l/n). 

The other off-diagonal it fields do not mix and stay with the mass parameter M 2 . Recalling the "88" component 
with mass M| 8 M 2 , one has n 2 — 2(n— 1) — 1 = (n— l) 2 pseudoscalar fields staying with mass M 2 , which corresponds 
to the breakdown of the diagonal symmetry U(n) to U(n — 1). Although the mass corrections of these "pions" are 
suppressed in the limit n — > 00, we write them explicitly, since they are essential for the renormalization of Eqs. |7]): 



Ml = M 2 - i 



^Y* = M 2 + -AM 2 . 
ri\ n n 



(12) 



The 2x2 coupled two-point matrix of the (K, Y 2 ) sector looks at n — 00 as follows: 



iG -i _ ( ^ - i^Y* - i ^v 8 



(13) 



Note that, in ( 13 1, the K index refers to a (x, j, k) type of generator, while the corresponding Y 2 carries the appropriate 
(y,j, k) type of index. There are also (n— 1) identical pairs with these indices reversed, where the off-diagonal elements 
of the 2x2 matrix receive an additional minus sign due to the totally antisymmetric nature of the f abc structure 
constants. The corresponding "kaon" mass is given by the zero of its determinant: 



M 2 K = M 2 +i x l]^Y 1 8 + ng 2 v : 
V in 



(14) 



The specific interest of the (K,Y 2 ) propagator matrix (13) is that it requires a nontrivial large n scaling for Yf and 
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us, if one wants to arrange finite contributions to Mf^ 



Xs 

I — ; 

\ n 



(15) 



It turns out that also the external source h% should scale as hg 
decouples from the EoS of vq and the SPEs of Yj and X: 



As a consequence, the vg condensate 



= nV2(h - M 2 v Q ), 



1 



T F (M 2 



92X 



(16) 



Here the notations y® = /n, x = X/n are introduced, and Tp(M 2 ) is the finite part of the tadpole integral of mass 
M: T F (M 2 ) = i J k (k 2 - M 2 )- 1 - - (M 2 - M%)Tf\ where M plays the role of the renormalization scale. The 

definitions of the quadratically divergent and the logarithmically divergent quantities were given in [16] . Their 
explicit form defines the renormalization scheme we used in the numerical part of our investigation. Here and below, 
we simply omit the divergent pieces from the tadpole integrals appearing in the EoS and SPEs. In the next section, 
where the renormalized effective potential is presented, the correctness of this formal procedure is demonstrated. 
When one makes use of the off-diagonal element of the 2x2 kaon-propagator G(k,y 2 ), calculated as the inverse of 



(13), one arrives, for the EoS (using the large n limit of the relevant structure constants) at 



= V2n (h 8 - M 2 xs + u 



Xs) ~ 9iXsTf{M 2 k ) 



(17) 



The corresponding saddle-point equation looks also quite simple, after one recognizes that ds a aTF(M 2 ) = by dg aa 
Therefore, one adds and subtracts from the contributions of the kaon-tadpoles the pion tadpoles and finds 



Vi 



2uoXs 



xl 



It f (m 2 k ) + It f (m 



(18) 



One can solve the equations of the 8 components, XSjUii with Vo,y®,x taken from their respective Eqs. (16), solved 



for a set of ho,g\, and <?2- These latter part of the calculation coincides with the solution presented in [IB], where only 
the vq condensate was considered. 



IV. THE RENORMALIZED EFFECTIVE POTENTIAL 

In case one finds nontrivial solutions for xs > yf also when h% = beyond the trivial xs — yf — solution, one has 
to decide which solution is preferred energetically, e.g., one has to evaluate the effective potential in the "light 7r" 
approximation described above. This means that, beyond the classical potential, written with the help of the auxiliary 
variables, also the contribution from the quantum fluctuations of the light particles, e.g., of the pions and kaons, should 
be included. The corresponding formal expression is the following, with 0(n) accuracy: 



v = v cl + v„ u 



Ye 



V 



quant 



quant ) 

M 2 v 2 + ^(x 2 + (y° 1 ) 2 )-2h Q v 



M Xs + 2(2/?) " 2h *Xs ~ W2g 2 ytxs{2v {) - X s) 



(n 2 - 2n) / ln(-p 2 + M 2 ) + 2n / ln(-p 2 + M 2 K ) 
v Jp 



(19) 



where V quan t is the standard expression of the one-loop part of the effective potential. It is important to note that, for 
the accurate eva luat ion of the quantum contribution of the pions, one has to make use of the 0(n) accurate value of 
M 2 , as given in (12). 

The divergences of V quan t can be eliminated by subtracting appropriate constants with divergences separated in 
terms of the auxiliary propagator i/(p 2 — M 2 ) in a way to have compatibility with the previously obtained finite Eqs. 
(16), (17) and (18). The counterterm part of the potential which absorbs these divergencies is the following: 



yc 



n 

Y 



(M 2 - m 2 )(TP - M 2 TD + ^(M 4 - m 4 )^' 



1 



(0) 



ln(-p 2 



M 2 ) 



-AM 2 



2 < /'-' i frr(2) /, /f 2 11^2^(0)% , 1 / 7. /r 2 i # 2\2/r(0) 



M 2 K - M ) (Tj + (M - M 2 )TD + \[M* K - M 2 ) 2 T^ 



(20) 



Although (20) is written in a rather compact form, it is a function of vq, xst m. 1 U\ 1 an< i x i an d is completely allowed by 
the formal requirements of renormalizability [T5]. Note also the important simplification AM 2 /2 + — M 2 = gix\- 
One can check explicitly that the derivatives of the renormalized potential (21 ) below, with respect to vq, xs, V\ , Vi and 




FIG. 1: Kaon/pion squared mass ratio as a function of the couplings ((71,32) for two different ho values. Kaons become lighter 
as ho increases, at the same time the region where nontrivial solutions are found narrows. 



x, reproduce the finite Eqs. (16), (17) and (18). One should also recognize that the second term of the right-hand side 
of the first line of (20) is field independent, and, therefore, gives no contribution to the equations for the one-point 
functions. It serves only for eliminating the overall divergence generated by the zero-point fluctuations of the pions. 

The quantum part of the renormalized effective potential (V quan t,R = V quan t + V ct ) in the renormalization scheme 
outlined above looks like 



Vq Uan t^R 



n 

n 
16tt 2 



{Ml 



M 4 ) + Mq(M 2 



2 (AA^-M 2 )Af 2 ln^ 



M 2 ) 
l 



2 M 2 



Mk In 



Ml 
M 2 



AM%M 2 + -(M 4 - M 



Ki 



(21) 



At this point, we have the complete renormalized effective potential, Vr = V c i + V quant ,R, and we can decide which 
solution is more favorable energetically Since we found two solutions for Xs, 2/i m a fixed background of vq, x, y values, 
we are only looking for the xs, V\ dependent parts of the effective potential to select the true ground state of the system. 
One can immediately see that only the C(n) part of Vr is what we have to study. It is also important to stress that, 
due to this fortunate effect, the selection of the ground state does not depend on the explicit value of n, since it appears 
as an overall multiplicative factor of every relevant term. In the forthcoming figures, we refer to V^ n ' as the 0(n) part 
of the full effective potential without the overall n factor. 

The fact that the 8 condensate contributes only to 0(n) means that the i>g contribution will be subleading relative 
to that of vq. Still, it belongs to the leading large n piece, since the next-to- leading-order contributions which were 
omitted from the very start would contribute O(n ). 



V. NUMERICAL RESULTS 



In [16], we mapped out the region of the three-dimensional parameter space ((71,(72,^0), where the scalar masses 
calculated from the solution of ( 16 ) turned out to be at least twice as large as the mass of the pseudo-Goldstone fields 
(pions). In the present investigation, we have chosen rather low values of ho (all dimensional quantities are measured 
in proportion to the appropriate powers of the absolute value of the renormalized mass). Note, however, that, in this 
paper, we have chosen Mq = — m 2 for the renormalization scale, different from the choice of [16]; therefore Fig. 2 of that 
paper can not be directly compared to this study. We just state that the results below are obtained in those regions of 
the parameter space, where the heavy scalar assumption holds (i.e., where the scalar masses are at least twice as large 
as the pion mass). This assumption will be checked if it also holds for the scalar/kaon mass ratio as well. 

Let us start the discussion with the mass splitting in the 7r sector - in particular, the kaon/pion mass ratio. The most 
detailed study was achieved for h s — 0, addressing the question of the v s (e.g., Xs) dependence of the effective potential 
after eliminating the yf dependence of Vr with the help of (18). We scanned through the allowed region at different 
values of ho- Our main result is that, in addition to the trivial solution v% — iyf — 0, also a nontrivial solution of (17) 
and ( 18 ) with positive squared masses of all tt fields was found in a large part of the region allowed by the heavy scalar 
mass assumption. In Figjlj the ratio M^/M 2 is displayed over the ((71,32) plane for 2 different values of h . One can 
see that, in a quite large region of the allowed parameter space, the mass ratio varies mildly, but the kaons are getting 
heavier when ho is lowered. Also, it can be observed that, for small g\ values, increasing gi induces very large kaon 
mass values, while, for larger g%, the same does not lead to significant changes of the mass ratio. 

The 0{n) part of the full renormalized effective potential Vr could be plotted in an extended interval around the 
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g,=5.5, g 2 =9.0, M =]m| 



V |n) /|m| 4 




Xs'M 



FIG. 2: The 0(n) part of the quantum effective potential divided by n as a function of the \s condensate. The real part of the 
potential is plotted with dotted lines for xs values where some negative mass squares lead to complex values of the potential. 
The curves from top to bottom follow the order of the ho values indicated in the upper-left-hand corner. 

nontrivial minimum, where all mass squares are positive. An example is shown in Fig(2j where one can see that the 
energy density of the nontrivial solution grows monotonically with ho increasing. [As already announced, the 0(n) 

part of the full potential divided by n is denoted as V'™', and the same holds for the classical part as \Q ■ .] There is 
some critical value /ig(<7i, (72), where the nontrivial local minimum disappears. For example, one can verify on Fig. [I] 
that, for the </i,<?2 values of Fig. [2j there does exist a nontrivial solution of the equations for ho = 0.005 but not for 
ho = 0.05. 

The dotted pieces of the curves represent the real part of the effective potential in regions where it becomes complex 
due to the fact that one of the mass squares becomes negative. This phenomenon is a well-known feature of the loop 
expansion [TTJ, [THl US] . Although the effective potential should be a real quantity by construction, in regions where 
the classical potential is chosen to be nonconvex, the loop expansion breaks down and produces an imaginary part. 
Mathematically, a nonconvex classical potential can be understood as an analytic continuation of a convex classical 
potential through the mass term: m — > im. The loop expansion creates a divergent asymptotic series, but there is no 
such theorem which would state that the analytically continued series corresponds term-by-term to the series of the 
analytically continued effective potential. In other words, the analytic continuation and the expansion in H are not 
interchangable operations. From a mathematical point of view, this is the origin of the appearance of the imaginary 
part. Nevertheless, the naive calculation leading to complex effective potentials has a physical meaning. First, it turns 
out that, if one wants to understand the effective potential at a <j) c point as the expectation value of the minimum 
energy density of the class of quantum states in which the field expectation value is 4> c (which property can be derived 
from the usual Legendre transform definition of the effective potential), then one has to take the real part of the 
loop expansion and apply Maxwell's construction. Furthermore, Weinberg and Wu showed [21] that even the naively 
calculated real part in itself contains physics: it is the minimum energy density of a class of quantum states with a C 
field expectation value, together with the more restrictive constraint that their wave functionals are concentrated on 
the uniform configuration <p c . The state chosen this way may differ from the state which minimizes the energy density 
of states with a 4> c field expectation value. In this case, the state in question is unstable, and the imaginary part of the 
potential can be interpreted as half the decay rate per unit volume. 

We have also plotted the curve of the classical part of the 0(n) piece of the potential in a typical point of the 
coupling space for two values of ho-, which demonstrates that, with the present renormalization scale Mo the quantum 
fluctuations only moderately modify the value of the potential; see Fig(3] In the whole allowed region of the coupling 
space, we found for the potential in the nontrivial minimum positive values; therefore, we conclude that, with h% = 0, 
there is no ho,vo,M 2 combination which would induce a spontaneously nonzero v% condensate into the true ground 
state of the system. The effect of the existence of such a metastable state might still show up in out-of-equilibrium 
situations. 

In a number of energetically "promising" coupling points, we have studied the effect of the application of a nonzero 
hg on the effective potential. In these cases, also the trivial solution shifts in proportion of the external field, and the 
effective potential can be displayed in a finite but small region around the shifted trivial minimum, too (i.e., no negative 
mass squares appear in its neighborhood) . If h is small, the value of the effective potential in the nontrivial minimum 
varies faster with h$ than near the shifted trivial minimum; therefore, one finds a critical external field hg(gi, (72) which 
drives the minimum characterized by a large v$ value to be the true ground state of the system. In Fig|4] we show the 
typical variation of the potential as a function of h$ for the value of (gi,g2, ho) chosen as in the left panel of Fig{3j In 
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h =0.005|m|3, g 1= 5.5, g 2 =9-0, M =|m| h =0.01|m|3, g 1= 5.5, g 2 =9-0, M =|m| 




0.1 0.2 0.3 0.4 0.1 0.2 0.3 0.4 

FIG. 3: The effect of the quantum fluctuations is only moderate, as indicated by comparing the Vj™' (dashed lines) with V 
The dotted part displays the real part of V^"\ where it becomes complex. 



h =0.005|m| 3 , g,=5.5, g 2 =9.0, M =|m| 




-0.005 1 1 1 1 1 1 1 

0.2 0.25 0.3 0.35 0.4 0.45 



FIG. 4: Smooth variation of V ^ with increasing hg displaying the symmetry breaking pattern Uv{n) — > Uv(n— 1). The curves 
from top to bottom follow the order of the applied h% values indicated in the upper-left-hand corner. 

Fig(5] we display the variation of the ratio M^/M^ as a function of h$ and g% for given g\ and ho values as taken in the 
nontrivial minimum. One can see that, apart from a jump to unity happening at lower g 2 values, the ratio is not varying 
too much. This ratio equals unity in the region where the nontrivial minimum is missing. Choosing higher values for 
ho involves smaller kaon masses (as it can also be seen in FigjTJ) , which points to the direction in the parameter space 
with a chance of phenomenological applications. One also has to note that, choosing other g\ values, the shape of the 
figure changes only mildly. The figure also shows that the mass ratio is increasing linearly with hg at fixed g 2 \ however, 
this change is almost negligible compared to its full magnitude, i.e., the relative change is very small. This means that 
the kaon mass is induced dominantly by the action of the vq condensate and not by the explicit symmetry breaking 
proportional to h&. This circumstance is rather different from what comes in the n — 3 flavor models from various 
perturbative discussions. 

Finally, we mention that it was checked that the scalar masses are also at least twice heavier than the kaons in 
the overwhelming part of the allowed parameter space, which means that our approximation based on the specific 
scalar/pseudoscalar mass hierarchy remains valid still in a wide range of parameters. 

VI. SUMMARY AND OUTLOOK 

Symmetry breaking patterns emerging in the ground state structure of the Ui,(n) x U^(n) symmetric matrix model 
were explored at leading order in the 1/n expansion. For the first time, ground states genuinely different from the broken 
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h =0.05|m| 3 , g 1= 2.0, M =|m| 

M K 2 /M, 2 




FIG. 5: The variation of M^/M^ with ho and <?2- The relative variation of the mass ratio is very small with respect to hg. Note 
that the abrupt jump downwards is occurring when the nontrivial minimum is missing. 

0(2n 2 ) symmetry were studied with a resummed perturbative approach, including the effect of the quantum fluctuations 
of the light (pseudo-Goldstone) fields. The effective potential of this approximation was consistently renormalized and 
used to compare energetically local minima corresponding to different symmetry breaking condensates. It was shown 
that in a large part of the coupling space, a condensate breaking Ul(ti) x Un(n) to LV(n) induces also the existence of a 
metastable local minimum with a further symmetry breaking Uv(n) — > Uv{n— 1). The effect of the explicit symmetry 
breaking sources on the local minima and the spectra was studied in detail. It was established that the application of 
an appropriate combination of the external sources drives quite often the Ul^ti) x Ur{ti) — > Uy(n — 1) breaking local 
minimum into the true ground state of the system. 

The different ground states are characterized by different pseudoscalar excitation spectra, also explicitly constructed. 
The mass ratio of the heavier modes ("kaons") and the remaining (n — l) 2 pseudo-Goldstone modes was found sensitive 
only to the external source ho proportional to the T° generator. The mild dependence on h^ ag T^ ag is rather different 
from what characterizes the kaon/pion ratio in perturbative treatments of the £/l(ti) x Un(n) — > Uv{n) symmetry 
breaking. If it turns out that one can build phenomenology on the nontrivial vacuum condensate, this observation 
would imply also rather different finite temperature behavior, compared to the perturbative descriptions. 

Although the transition in the vq condensate in itself would be much similar to that of the 0(2n 2 ) symmetric model, 
the insensitivity of the location of the nontrivial minimum in v^ 9 with respect to h„ a9 (see Fig[5]) hints to a first-order 
transition. We observed that it is not h°^ a9 but rather v which determines the value of vf™ 9 . The tightly connected 
variation of the two condensates probably leads then to a phase transition with both fields going through a discontinuous 
change of similar amplitude. This expectation complies well with renormalization group considerations which classify 
the restoration of the broken symmetry in this model quite differently than in models displaying symmetries of the 
orthogonal group H3] . An obvious next step is the study of the finite temperature symmetry restoration from this 
nontrivial ground state. 

Further consolidation of the discovered structure requires the investigation of the fluctuation effect of the heavy 
scalar fields. On a general basis, one expects the modification of the different coefficients in different EoS, SPEs, and 
self-energies by terms of 0(M 2 /M 2 ) (where M 2 is one of the scalar masses), which would just deform slightly the 
coupling region where the nontrivial vacuum structure is self-consistently present. 

Our final goal is to make use of this solution at n = 3, when also the contributions from the effective term reflecting 
the axial anomaly (the 't Hooft determinant) would be included perturbatively using the large-n propagators of the 
different fields. The comparison with the treatments of the n = 3 three flavored meson model in different perturbative 
calculations ^4-28] and functional approaches, which go beyond perturbation theory partly based on the exact 

renormalization group techniques |29H31j . could give some insight into the importance of strong coupling effects. The 
further exploration of the proposed ground state will proceed along these three lines. 
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